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Abstract 



We test the capability of state-of-the-art lattice techniques for a precise deter- 
mination of relativistic corrections to the static interquark potential, by use 
of SU(2) gauge theory. Emphasis is put on the short range structure of the 
spin dependent potentials, with lattice resolution a ranging from a ~ 0.04 fm 



(at (3 = 2.74) down to a « 0.02 fm (at (3 = 2.96) on volumes of 32 4 and 
48 4 lattice sites. We find a new short range Coulomb-like contribution to the 
spin-orbit potential V[. 
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I. INTRODUCTION 



Quarkonia spectroscopy provides a wealth of information and thus constitutes an im- 
portant observational window to the phenomenology of confining quark interactions. It has 
been known for a long time that purely phenomenological or QCD inspired potential models 
offer a suitable heuristic framework to understand the empirical charmonium (J/i/j) an d 
bottomonium (T) spectra 

On a more fundamental level, one would prefer to start out from the basic QCD La- 
grangian to solve the heavy quarkonia bound state problem. Two alternative strategies lend 
themselves for reaching this goal: (i) direct extraction of the bound states on the lattice 
from an effective nonrelativistic lattice Lagrangian approximation (NRQCD) ||, (ii) use of 
an effective nonrelativistic Hamiltonian framework through the intermediary of potentials 
determined from lattice QCD. 

Considerable efforts have been made recently to determine the quarkonia spectra within 
the NRQCD approximation of QCD |J; the notorious technical problems to determine ex- 
cited states in the Euclidean formulation have been tackled with remarkable success. In 
the alternative Schrodinger-Pauli setting, the technical problems are shifted towards the lat- 
tice determination of 1/m 2 corrections to the potential. The spin dependent (sd) ||[7] and 
velocity dependent (vd) [§] contributions need to be extracted from (Euclidian) time asymp- 
totia of rather complex observables that require renormalization and must obey constraints 
following from Lorentz symmetry ||[5|]. 

First attempts to compute the relativistic corrections to the static potential on the lattice 
have been pioneered in the mid eighties flOHTBl. In the meantime tremendous improvements 
have been achieved both in computational power and methods. The central potential has 
been determined with high accuracy in quenched QCD |l4|-p^| and, more recently, in full 
QCD with two dynamical flavors of light Wilson sea quarks JH]]. In view of the general 
interest in the potential formulation of the meson binding problem, renewed effort should 
be made to unravel the structure of sd and vd potentials. This will provide us with a 
better understanding of the structure of the interaction in the intermediate distance regime 
0.15 fm < r < 1 fm which is of tantamount importance to the binding problem. 

As a first step within this program we shall present in this paper a high statistics study 
of the spin dependent forces in SU(2) gauge theory. Though the two color formulation 
will not yet allow to proceed to spectrum calculations we would expect the key features of 
gluodynamical confinement to be revealed. In a follow-up paper JE^] (referred to as II) we 
shall apply our techniques to the SU(3) case. 

The present article is organized as follows: in Section |J, we provide an introduction into 
the Hamiltonian formulation of QCD binding problems, the expected theoretical scenario of 
potentials as well as the lattice observables from which to determine them. In Section fH\ , 
useful lattice techniques will be elaborated. In particular we shall discuss the systematic 
uncertainties of the approach. The resulting SU(2) gauge potentials will be presented in 
Section [TV]. 

II. THE HEAVY QUARK POTENTIAL 
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A. Hamiltonian formulation of the meson binding problem 



Starting from a Foldy-Wouthuysen transformation of the Euclidean quark propagator in 
an external gauge field, the asymptotic (T — > oo) expression (W(R, T)) oc exp(— Vo(R)T) for 
expectation values of Wilson loops can be derived. V (r) = a _1 Vo(i?) denotes the potential 
between static quarks, separated by a distance r = aR. R and T are the spatial and 
temporal extents of the (rectangular) Wilson loop. By perturbing the propagator in terms 
of the inverse quark masses m^ -1 and m^ 1 around its static solution, one arrives at the 
semi-relativistic Hamiltonian (in the CM system, i.e. p = Pi = — P2 and L = L x = L 2 ), 




2 4 
p p 



+ V {r) + V ad {r, L, S 1} S 2 ) + V vd (r, p), (1) 



2rrii 8m'} 

where the potential consists of a central part, sd ||[7j and vd || corrections^, 

LSi + LS^ \ V>(r) + 2V;(r) 
m\ m\ ) 2r 




with 



and 



% = rj p£-*i ( 3 ) 



K d (,P)4(^ + ^)(vV oW + VV a ( 



1 , „ , A 1 



{Pi,Pj, Sij } w j + ^2 2 {pj , Pj , Ty } w j 

mim 2 J fe=1 m fc 



(4) 



with S^- = 5„y 6 (r) - RijV c (r) and 7y = (^V^r) - R^V^r). The symbol {-, -, -} Wey i = 
H''{'''}} denotes Weyl ordering of the three arguments. V{ - V4 are related to spin-orbit 
and spin-spin interactions, Vb-V e to orbit-orbit interactions and the Darwin-like term that 
incorporates V 2 V a modifies the central potential^. V[ - V4 and V 2 V a - V e can be computed 



1 Here, we just state the classical values for the couplings to the potentials. Radiative corrections 
give rise to small deviations from these tree level results fpTOfl - The corresponding one loop coeffi- 
cients have recently been determined in the framework of heavy quark effective theory (HQET) in 
Ref. |2C]. In the unequal mass case additional contributions appear, whose tree level coefficients 
are zero, that can also be parametrized in terms of V[ an V^- 

2 In following the convention of Ref. || we have included this term into V vd although it does not 
explicitly depend on the velocity. 
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from lattice correlation functions in Euclidean time of Wilson loop like operators. Pairs of 
the potentials are related by Lorentz invariance to the central potential ||J|: 



VZ(r)-V{{r)=V£{r), (5) 
Vb {r) + 2V d {r) = r -V'{r) - l -V Q {r), (6) 

Vc (r) + 2V e (r) = - r -V^r), (7) 

such that only six sd and vd potentials turn out to be truly independent. 

In the present SU(2) investigation, we restrict ourselves to the spin dependent terms. 
Note, that the Hamiltonian Eq. ([!]) contains dimension six operators and thus is not renor- 
malizable. For this reason the theory — being truncated at order 1/m 2 — is only an effective 
one with validity range of small gluon momenta (compared to the heavy quark masses). This 
very fact gave rise to a discrepancy between the Eichten-Feinberg-Gromes formulae Eq. (^) 
and perturbative expansions pTj|l9| in powers of the coupling, g, where additional logarith- 
mic mass dependencies occurred from dimensional regularization. The underlying problem 
is now solved and one loop matching coefficients between the effective Hamiltonian and QCD 
have been obtained from HQET [pO . 



B. Expectations on sd potentials 

In addition to the exact constraint, Eq. (|5|), derived by Gromes |J, some approximate 
relations between the sd potentials are anticipated from exchange symmetry arguments. We 
start from the standard assumption that the origin of the central potential is due to vector 
and scalar-like gluon exchange contributions. Given the fact that a vector-like exchange 
can at most grow logarithmically with r |[22|| , the nature of the linear part of the confining 
potential can only be scalar. As we will see, V^r) is short ranged, such that the confining part 
only contributes to V[(r). This leads us to expect V^'(r) to be purely vector-like. Under the 
additional assumptions that pseudoscalar contributions can be neglected and that V[ does 
not contain a vector-like contributions, one ends up with the scenario of interrelations ||, 

V 3 (r) = ^-V!(r), (8) 
Vt(r) = 2VV 2 (r), (9) 

which of course has to be in agreement with leading order perturbation theory. However, 
Eqs. ©— © hold true for any effective gluon propagator that transforms like a Lorentz 
vector. In principle, to this order in m _1 , V3 and V4 could contain additional pseudoscalar 
pieces. 

Tree-level continuum perturbation theory (see Appendix) yields the following expecta- 
tions for the central, spin-orbit and spin-spin potentials: 

V {r) = -p (10) 
V{{r) = 0, (11) 
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Vlir) = (12) 

Vs(r) = ^ (13) 

V 4 (r) = 8vre5 3 (r), (14) 

where e = Cfo> s and Cp = 3/4 for SU(2). Note, that combining the perturbative result on 
V^'(r) with the Gromes relation, Eq. (|5|), and a funnel type parametrization of the central 
potential [0, 

V (r) = -- + Kr, (15) 

the first spin-orbit potential should take the form 

V{{r) = -k. (16) 

In the Appendix we shall also derive the corresponding tree level results for the lattice 
potentials. We will show that an exact lattice relation between V{ and V 2 ' and a linear 
combination of the central potential, taken at different lattice coordinates, in analogy to the 
Gromes relation does not exist. However, the relation should be retrieved in the continuum 
limit and approximately hold within the scaling region on the lattice. 



C. How to compute the potentials 

In the potential picture, the quarks interact instantaneously through a potential that 
only depends on the distance, spins, and velocities of the sources, Eqs. ([I]), (0) and (|]). 
All time dependence has been separated and implicitly included into coefficient functions of 
various interaction terms, the so called sd and vd potentials. These can be computed by a 
non-perturbative integration over gluonic interactions. Therefore, the sd and vd potentials 
incorporate a summation over all possible interaction times, t. One obtains the following 
expressions in terms of expectation values in presence of a gauge field background for the 
sd potentials 

~R, 

§V£(R) = e ijk lim rtZtt((S f (0,0)4(R,t)» w , (18) 
RijV^R) = 2 lim fdt [(( J B i (0,0)4(R,t))) w (19) 

T ►OO Jq L 



-V((R) = 2e ijfc hm f dtt ((5,(0, 0)^(0, t))) w , (17) 

J 



^<<B(0,0)B(R,t)>> 



w 



V 4 (R) = 2 Jim J^dt «B(0,0)B(R,t))) w . (20) 



We have put the expressions into a form that is more suitable for lattice simulations. Via spectral 
decompositions of the underlying correlation functions, equality of the expressions to those of 
Refs. |(|[?]] can easily be shown. 
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With a — > 0, the above potentials should approach their continuum counterparts and ro- 
tational invariance should be restored, V{(R) = a 2 V{(aR), V^'(R) = cPV^aR), V 3 (R) = 
a?V 3 (aR) and K 4 (R) = a 3 V 4 {aR). 

The operator ((FiF 2 ))w is defined as follows: 



«^2» 



(Tr V [exp (ig J dw dx^ A„) 



F 1 F 2 \) 



w 



IdW A")]) 



(21) 



(Tr V [exp (ig f 9 

where dW denotes a closed path (the contour of a Wilson loop W(R,T)). V denotes path 



ordering of the arguments. Following Ref. |T3[ we have chosen the following discretized 
version of Eq. (pTF) for the case of two color field insertions: 



(F l F 2 ) 



w 



(V 


w (nx-nl)^ 


-ut)_ 


)(W) 


(V 


w (ri! + u{) 


)(v 


w (n 2 + r4)~ 


> 



(22) 



where the subscripts 1 and 2 represent the multi-indices (ni,/^^) and (n 2 , fi 2 , v 2 \ re- 
spectively, rti are integer valued four-vectors. The above ratio of "eared" Wilson loops is 
visualized in Fig. |l| In order to avoid imaginary phases and factors g 2 a A from Eqs. (|17|)-(p0l), 
we use the following conventions for the electric and magnetic fields: 



hp 



ga 2 F, 



hi' 



E; 



B; 



c^-ijkFjk- 



(23) 



In what follows, we have chosen IL^n) to be the average of the four plaquettes, enclosing 
the lattice point n, 



n^(n) 



[P„, v (n) + P^(n) + P- M ,-„(n) + P-w( n )] 



(24) 



with 



Pw(n) = U„{n)U v (n + p)lA(n + v)Ul(n) 



fi). 



(25) 



This choice of II makes Eq. fl22| ) correct up to order a 2 , the discretization error of the Wilson 
action, used for generating the gauge field background. 

In practical computation, the temporal extent T of the Wilson loop W within Eqs. (|l~7|) - 
(pOj) is adapted according to the formula T = t + 2At. We choose to keep the minimal 
distance, At, between the "ears" and the spatial closures of the Wilson loop fixed. Note, 
that, strictly speaking, Eqs. (p!7|)-(|20|) apply for the limit At — > oo. Our coordinates are 
such that the Wilson loop VT(R, t + 2 At) extends from n 4 = —At to n 4 = t + At into 
the temporal direction and from n = to n = R into the spatial directions. Physically 
speaking, At represents the time we allow the gluon field to decay into the ground state, 
after (before) creation (annihilation) of the gg-state. Hence this deexcitation time must be 
considered as an important control parameter of our measurements. 
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III. LATTICE SIMULATIONS 



A. Simulation parameters 

Since, apart from V{, all sd potentials are expected to be short ranged, we will focus our 
attention onto short distance properties, i.e. aim at the smallest lattice resolution possible on 
our computers. The present simulations have been performed on V = L^L T = 16 4 , 32 4 and 
48 4 lattices at (3 = 2.74 and (3 = 2.96 which correspond to lattice spacings a ~ 0.041 fm and 
a ~ 0.022 fm, respectively (Table [I]). The number of independent Monte Carlo configurations 
^conf, generated at each set of parameters, is included into the Table. The above physical 
scales have been adjusted such that the string tension comes out to befj y/n = 440 MeV. 
Using such small physical volumes, finite size effects (FSEs) have to be investigated. This 
will be done by comparing results obtained on a 16 4 lattice with 32 4 results at (3 = 2.74. 
Finite lattice resolution (i.e. finite a) effects are investigated by relating results, obtained at 
the two different values of the coupling. 



B. Updating algorithm 

The numerical calculations are performed on lattices with hypercubic geometry and 
periodic boundary conditions in all four directions. Throughout the simulation the standard 
Wilson action 

Sw = -P E ^TrP^(n) (26) 

n,iX>f 

with (3 = 4/g 2 has been used. 

For the updating of the gauge fields, a hybrid of heatbath and overrelaxation algorithms 
has been implemented p3 |. The Fabricius-Haan heatbath sweeps |24[] have been randomly 



mixed with the overrelaxation step with probability 1/14. The links have been visited in 
lexicographical ordering within hypercubes of 2 4 lattice sites, i.e. within each such hypercube, 
first all links pointing into direction 1 are visited site by site, then all links in direction 2 etc.. 
After at least 2000 heatbath thermalization sweeps, measurements are taken every 100 or 
200 sweeps at (3 = 2.74 and (3 = 2.96, respectively, to ensure decorrelation. We find no signs 
of any autocorrelation effects between successive configurations within any of the measured 
observables. 



C. Link integration 

Statistical fluctuations have been reduced by "integrating out" temporal links that ap- 
pear within the Wilson loops and the electric ears analytically, wherever possible. By "link 
integration" we mean the following substitution [[251: 



4 This scale can only provide a rough orientation as here we are simulating SU(2) gauge theory, 
not full QCD. 



with 



and 



U*(n)^W*(n) = )^ dUe _ pSnAu) (27) 



SnAU) = "^Tr (UFfrnj) , (28) 



F» = J] U v (n)U^n + *>)E#(n + A). (29) 

Wi(n) is in general not an SU(2) element anymore. 

In this way, time-like links are replaced by the mean value they take in the neighborhood 
of the enclosing staples F±(n). Only those links that do not share a common plaquette can 
be integrated independently. 

In case of SU(2) gauge theory, W±{n) can be calculated analytically, 

w ' {n) = (30) 



where f^(n) = ydet(F M (n)). /„ denote the modified Bessel functions. 

D. Smearing 

In order to achieve a satisfying overlap between the quark- ant iquark ground state and 
the state created by the spatial parts of the Wilson loop, we have applied a smearing pro- 
cedure [11,^31 by iteratively replacing each spatial link Ui{n) within the Wilson loop, by a 



'fat" link, 

UAn 



{n) ^N^a Uiin) + £ U^U^n + j)U]{n + i) j , (31) 

with the appropriate normalization N and free parameter a. We find satisfactory ground 
state enhancement with the parameter choice n iter = 150 and a = 2. In Fig. |2| the resulting 
central interquark potentials are displayed. In Fig. |3], the corresponding ground state over- 
laps Co(R) are shown as a function of the source separation at (3 = 2.74. As can be seen, 
all overlaps are well above 0.9. 



E. Spectral decomposition 

In this section we will discuss the control of excited state contributions at finite deexci- 
tation time At. As explained above, the spatial transporters within the Wilson loops have 
been smeared to suppress such pollutions from the very beginning, allowing to work with 
moderate values of At. This is vitally important, as statistical errors increase with the size 
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of the Wilson loop. In practice, we have decided to keep At fixed at a value suited for 
sufficient deexcitation and to increase the temporal extent of the Wilson loop according to 
T = t + 2At, i.e. with the separation t between the two ears. We found At = 2 to be 
appropriate. 



Previous authors ||lT| , |r2f have replaced the integrals over interaction times by discrete 
sums. This results in cut-off errors due to the finiteness of r as well as additional order 
a 2 integration errors. Both sources of systematic uncertainties can be significantly reduced 
by exploiting transfer matrix techniques. For illustration of such techniques we start from 
the static potential which can be computed from Wilson loops: at Euclidean time t — 0, a 
creation operator, 

rl = q(0)U(0 - R)g f (R) (32) 

with a gauge covariant transporter {7(0 — > R) is applied to the vacuum state |0). g*(R) 
creates a heavy quark spinor at position R. The qq pair is then propagated to t = T by 
static Wilson lines in presence of the gauge field background, and finally annihilated by Tr. 
A spectral decomposition of the Wilson loop exhibits the following behavior (T = e~ aH 
denotes the transfer matrix, T\n) = e~ aEn \n)): 



P W = (W(R,T)) 



Tr (T R T T TlT L T- 



Tr (T Lt ) 

= \ T EK^|rR|^R)| 2 e-^ (R)T e^^- r ) (33) 

Z^m c m,n 

= ]T K(R)|V^ R ) T x (l +0 (e-^- T >)) 

n 

with 

d n (R) = (0|r R |n,R> , C n (R) = K(R)| 2 - (34) 

\n, R) is the nth eigenstate of a qq pair, separated by a distance R. |<i n (R)| 2 takes a 
positive value, whenever such a state has a non-vanishing overlap to the creation operator 
T R , applied to the vacuum. |n) is the nth pure glue eigenstate (glueball). V^(R) denotes 
the nth excitation of the qq potential and the vacuum energy Eq has been set to zero. E\ 
is the mass gap, i.e. the mass of the Af glueball. Due to its heavy mass and L T ^> T, such 
back-propagating terms can be neglected to high accuracy. The overlaps are normalized, 
J2 n Cn(R) = 1) such that a ground state overlap Co(R) ~ 1 implies C n (R) 1 for all n > 0. 
Co(R) can be increased by optimizing the path combination {7(0 — > R) within the creation 
operator Eq. (|32|) (smearing). 



Now, let us define the operators, 

^ = gl(ni-nt) g t , ft = g i(n i + nt) 9 t > (35) 

i.e. qlliq^ = Gi + for the chromo field insertions within the nominator and denominator 



of Eq. (p2"D. i = (rij,/ij,z/j) denotes the position of the insertion as well as the color field 



component. In the present case, rij either takes the position or R. Let 
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r mn (R) = (m, R\^\n, R) , <4 n (R) = (m, R|ft|n, R). 



(36) 



Hermiticity of ^ and Q { implies = f nm and ^ n = 

A spectral decomposition of a Wilson loop, W(H,T) with ears T\ and JF 2 inserted 
at times n 4 = and n 4 = t, respectively (The Wilson loop extends from n 4 = —At to 
n 4 = t + At.), yields (where we have neglected the back-propagation of pure glue states): 



D 12 
W 



V 



w(r,t) (iix-ni) (n 2 -n 2 



Tr (rT^^Tt^T^T^-T 



Tr (T Lt 



(37) 



(AVi+AV n )At e -AV m t 



l.m.n 



with AV n — V n — Vq. For better readability, the R dependence of Tr, V^(R), d n (R) as well 
as /^ n (R) has been omited from the above equation. 

Next, let us investigate the behavior of the Wilson loops with one ear, appearing in the 
denominator of Eq. (^2|). For an insertion of Qi at time n 4 = or n 4 = t we find, 



Q 



w 



t 

i 

L T -T 



w(r,t) (ir + n 
Tr (rT A *^T t+A 'rtT 

Tr (T Lt ) 



(38) 



In combining the above expressions we finally obtain 

Pw^w 



w 



QwQ 



2 

W 



^D^e-^ x 1 + E^e~ Av ^ At 



(39) 



with 



and 



D 



12 



R- e (fomfmo) 

PooPoo 



E 



12 



Re [(dx/do) (fLflo + fLfD\ ^ [{d 1 /d )g\ \ Re [(di/d )^i ] 



R- e (/om/mo) 



Poo 



^00 



(40) 



(41) 



Again, all constants are understood to depend on R. 

Note, that unwanted excited state contributions are suppressed by the ratio di/do as well 
as by e~ AVlAt . The smallest value of At that appears within an integral over interaction 
times will determine the reliability of the result. The bosonic string picture expectation 
is AVi(R) = tt/R ||28|| . This has been qualitatively confirmed in numerical studies p9| , |30 
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such that At = 2 yields an excited state suppression by a factor exp(— 4n/R) for large R, 
in addition to the ratio di/d . 

The creation operator projects only onto states within the A\ g representation of the 
appropriate symmetry group [p9[| . The lowest continuum angular momentum to which 



it couples is L = 0. The hybrid (L = 1) state E u is the next excitation. The combination of 
magnetic ears we use, applied to a qq state within the Ai g representation, results in a pure 
E u state which has no overlap with A\ g , such that Dq 2 = within all correlation functions 
of interest, i.e. all correlators decay exponentially with Euclidean time t. This does not hold 
true for some of the vd potentials. 

Unlike in the case of the central potential, where the sum of the (non- negative) overlap 
coefficients C m is normalized to one, the are not normalized and can be negative. 
However, due to invariance under time inversion, J2 m ^i? = in case of V[ and V{ since the 
correlation function has to vanish at t = 0. In combining Eq. (|39D with Eqs. ( |19D or (p0|), 
we obtain for V 3 or V4: 

POO A 

V 3A oc E / dtD^e"^ =E^- (42) 

m>0 m>0 

with appropriate color field positions ni = 0, n 2 = R and components /ii, U\, ^2, ^2- Eqs. ( |17D 
and QT|) yield: 

K, 2 oc E rdttD^e-^ = E t^V- (43) 

m>0^° m>0 f AV m ) 



From the above formulae it is evident that not only eigenvalues of the transfer matrix 
but also amplitudes enter in the computation of the sd potentials, which illustrates where 
renormalization constants come in, which are effectively removed by the multiplication with 
(S'ooS'oo) 1 [Eq- (pEPPl- This factor originates from the denominator of Eq. (|22). The pa- 



rameters and AV m can be fixed from a fit to Eq. (|55|). The hybrid potentials can 
in principle also be determined independently [p9[] . We leave this for future high precision 
studies on anisotropic lattices. For the time being, we evaluate the integrals Eqs. (|T7| ) -(|20| ) 
numerically where our interpolation method has been inspired by the multi-exponential 
result of the spectral decomposition, Eq. (|39|) . 



F. Interpolation procedure 

In order to determine the sd potentials, one must evaluate integrals over correlation 
functions [see Eqs. (P~7|) — (POD ] that depend on the interaction time t in a multi-exponential 
way, Eq. ([39|). In the following, C;(t) will denote the two point function which has to be 
integrated out in order to determine at a given value of R. For i = 1,2, Cj(t) will be 
weighted by an additional factor t [Eqs. ([T7|)-(P^D1. Two different methods of interpolating 
Ci(t) in between the discrete t- values have been adopted: 

1. Wherever the quality of the signals allowed for reasonable fits, the data has been fitted 
to a two-exponential ansatz, 
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C{t) = D x exp(-AV"it) + D 2 exp(-AV 2 t). 



(44) 



This amounts to three parameter fits in case of C\ and C 2 (Di + D 2 = 0) and four parameter 
fits for C3 and C4. Unfortunately, stable fits have only been possible for the region of small 
R where we do not necessarily expect the hybrid potentials to agree with string model 
predictions. 

2. Alternatively we have performed a local exponential interpolation: 

C i (t')=C i (t)exp[-B i (t)(t'-t)] , = ln (c^TT)) ( 45 ) 

for t < t' < t + 1 and Ci(t)Ci(t + 1) > 0. Due to the multi-exponential character of 
the correlation function (or statistical fluctuations) the sign might change within the given 
interval. Thus, for Ci(t)Ci(t + 1) < 0, we interpolate linearly, 

C t (t') = d(t) + [Ci{t + 1) - Ci(t)](f - t). (46) 

For C\{t) and C 2 (t) quadratical interpolations are performed within the interval < t' < 1 
to account for Ci(0) = C 2 (0) = 0, where we demand continuity of the interpolating function 
and its derivative at t — 1. 

The interpolation procedure is illustrated in Figs. ^-|7| for some examples. As can be 
seen from Fig. [| where C\(t) is displayed for R = 10, the quadratic interpolation for t < 1 
might differ significantly from the fit. However, this region does only weakly contribute to 
the potential since C\(t) as well as C 2 (t) are weighted by an additional factor t. In Figs. |5| 
and |, C 2 (t) is displayed for R = 2 and R = 4, respectively. C±{t) at R = 2 changes its sign 
as can be seen from Fig. [?[ 

All statistical errors have been bootstrapped. For each potential V^'\ numerical integra- 
tion has been performed up to a value t = Tj, with Tj chosen such that the result is stable 
(within statistical accuracy) under the replacement T{ — ► T{ — 1 for all R. Subsequently, 
systematic cut-off errors have been estimated from the exponential tail of fits to large t 
data points but came out to be negligible in all cases when compared to the statistical er- 
ror from the numerical integration up to Tj. Whenever the deviations between the fit and 
interpolation results turned out to be significant (see e.g. Fig. |7|), we have included them as 
a systematic uncertainty into the final error on the potential value. 



G. Renormalization and matching 

The sd and vd potentials are computed from amplitudes of correlation functions rather 
than from eigenvalues of the transfer matrix. This gives rise to renormalizations in respect 
to the corresponding continuum potentials. A different way to illustrate the necessity of 
renormalization is the fact that the color electric (magnetic) ears, 

F, v = ga 2 F, v = 1 (u^ - njj (l + 0(a 2 )) , (47) 
explicitly depend on the lattice scale a and, therefore, discretization. 
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As in the low energy regime of interest the renormalization constants are likely to re- 
ceive relevant high order corrections, we apply the non-perturbative HM renormalization 
prescription |U| [cf. Eq. fl22|) and Fig. [I[ that is geared to "dividing out" the two amplitudes 
(?q and #q [Eqs. (fhj), AUD I from the naive lattice two point function. Terms that do not 
depend on a dimensionful parameter (which is the distance between the sources, R, in the 
case of interest) will give rise to divergencies in the continuum limit. The idea behind the 
HM procedure is to systematically remove such terms from the correlation functions in order 
to arrive at residual renormalization constants that are close to one. 

In terms of perturbation theory, one can classify the related diagrams into pure self in- 
teractions within the ears, pure interactions within the Wilson loops, interactions between 
ears and the loop and — to higher orders — mixes thereof. By our renormalization proce- 
dure all divergent diagrams that consist only of ear or Wilson loop self interactions, cancel 
identically. In addition, many more complicated contributions are removed. Within V% - 
V4, all divergencies of orders g 6 and less vanish while in case of V[ this holds only true up 
to orderQ g 4 , such that the renormalization constants involved will only differ from identity 
on a three loop [1 + 0(g 6 )] or two loop [1 + 0(g 4 )} level, respectively. 

The spirit of the procedure is close to the one of Lepage-Mackenzie who suggest to 
construct "tadpole improved" operators |JT[]. In fact, one finds the HM renormalization 



and the tadpole improvement prescriptions to render near-equal results (within 2 %, as 
compared to overall renormalization effects of 60 %, at our present /5-values). The remaining 
difference between the "tadpole improved" operator and the HM renormalized counterpart 
can be explained in terms of more complicated diagrams involving interactions between the 
ears and the time-like parts of the Wilson loop. We take the small size of this difference as 
an indication that other similar higher order terms, which we have not been able to cancel 
out completely, can be neglected. 

Direct numerical checks of the accuracy of the approach are possible in two ways, namely 
(i) by varying the lattice resolution a and testing scaling of the results^] and (ii) by comparing 
the data with the prediction on V 2 ' — V{ from the Gromes relation, Eq. (|5p, between spin-orbit 
potentials and the central potential (which does not undergo renormalization), 

V^r) = Vl rcn (r)-Vl rcn (r). (48) 

In Fig. [8] we check our data on V 2 ' — V[ in units of the string tension k against the force, 
obtained from a fit to the central potential, Vo(r) [Eq. (|5"1"D, Table |TJ. As can be seen, the 
two data sets scale beautifully onto each other and reproduce the central force up to lattice 
artefacts at small r. 



5 Within V{, gluon exchanges between the two ears that are not cancelled by the denominator of 
Eq. ( |22"|) contribute to the self energy. 

6 However, due to the running of the matching constants to the full theory with the lattice scale, 
residual scaling violations of about 2.7 % for V{, and V% and 5.5 % for V3 and V4 are expected after 
having rescaled the potentials by relative factors as large as 02.74/°!. 96 ~ 3-2 and ai.74/ a 2.96 ~ ^-7, 
respectively. Needless to say, that we cannot resolve such tiny effects from our lattice data. 
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Renormalization is not a pure lattice problem in this case. By truncating the 1/m 
expansion of the SU(2) Lagrangian at a given order, the ultra-violet behavior is altered in 
respect to the full theory. Therefore, the resulting effective Lagrangian has to be matched to 
full SU(2) gauge theory at a renormalization scale \i = ir/a giving rise to coefficients m), 
connecting an SU(2) potential at the heavy quark mass m, Vi(r;m) to the corresponding 
potential, computed in the framework of the effective theory at scale /i, e.g. V^(r;m) = 
Ci(n, m)Vi(r; /i). This problem, which becomes visible beyond the tree- level, has been treated 
in a systematic manner for sd potentials by Chen et al. [f2(J. The Gromes relation still 
remains valid [V '(r) = V% (r; jj) — V{(r; fi)], and constrains the matching coefficients whose 
one loop values have been derived in the Reference. Unlike the renormalization constants 
that relate the lattice potentials to the continuum counterparts, in this case, we can rely on 
perturbative results, with the one loop correction being an effect of just a few per cent. 



H. Finite volume effects 

It is well known from previous simulations of SU(2) gauge theory (cf. Ref. ]2E|] ) that 
finite size effects (FSEs) on the central potential, computed from Wilson loops, and the 
string tension are negligible within the numerical accuracy of typical lattice studies, as long 
as the spatial lattice extent is kept larger than L a « 1 fm (and L T > L a ). However, since the 
critical temperature of the deconfinement phase transition corresponds to a scale of about 
0.65 fm P2]| , the slope of the potential decreases rapidly as the extent of the box is reduced 
below about 0.8 fm. Nonetheless, one might expect the sd potentials V{ - V4 to be affected 
much less by the infrared cut-off than the central potential, as they are of much shorter 
range than the latter [see Eqs. (0), (0)-(0) or Eqs. (|M7|) -([K40|) of the Appendix]. 

All results at (3 — 2.74 have been obtained on a volume of 32 4 lattice sites which is 
comfortably large. In addition, we have performed simulations on a 16 4 lattice, whose 
physical extent approximately corresponds to that of a 32 4 lattice at (5 = 2.96 (Table |), 
to estimate FSEs on the individual potentials. At (5 = 2.96, we have determined the string 
tension and parametrization of the central potential on a 48 4 lattice which is sufficiently 
large to avoid serious FSEs. The extracted scale a = 0.022 fm is in agreement with the 
expectation we have from the data collected in Ref. [po ]. As expected, we find significant 
FSEs on the central potential obtained on a 32 4 lattice at this /3-value as can be seen from 
Fig. [| The slope is reduced by about 35 % on the small volume while the string tension on 
a 16 4 lattice at /3 — 2.74 comes out to be smaller by as much as 50 %, compared to V = 32 4 . 

At j3 = 2.96, the computer memory forced us to restrict ourselves to the smaller volume 
(V = 32 4 ) for the sd potentials. In the following, we compare results obtained on the 16 4 
lattice to V = 32 4 results at (3 — 2.74 to justify the assumption that the j3 = 2.96 sd 
potentials are not seriously affected by FSEs. As can be seen from Figs. |10|-13, for all sd 



potentials FSEs are statistically insignificant. Differences between the constant long range 
contributions to V[, which we do expect, are hidden within the large statistical errors of the 
{{Bi(0, 0)Ej(0, t)))w correlators. We conclude that though a spatial lattice extent of 0.7 fm 
is too small to extract the infinite volume central potential, it suffices for extracting V 2 ' - V4 
as well as the short distance contribution to V[. 
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IV. PHYSICS RESULTS 



A. Getting started: the central potential 

The lattice potential data Vo(R) has been computed from smeared Wilson loops by the 
method described in Ref. f26|. Our general strategy is to derive continuous parametrizations 
of the lattice data points which will enable us to construct derivatives and compare the results 
to theoretical expectations, such as the exact Gromes relation [Eq. (|5|)] or the approximate 
relations, Eqs. (|8])-(|9]). The continuum tree level results on the central and sd potentials 
have been presented in Section |TTB| [Eqs. (|T(i|)- (lT4]) l and the corresponding lattice expressions 
are derived in the Appendix [Eqs. (|A§), ( [MTD , ( |A25|) -(|A2^) and (fA2D|)], such that we can 



correct the lattice data for the differences before attempting to fit them to a continuous 
parametrization. 
Let 

Vcon(R) = Vo(R) - gSV (R). (49) 

with 

5V (R) = ±--4tcG l (R) (50) 
K 

be the tree level corrected central potential. The lattice gluon propagator Gl{R) is defined in 
the Appendix [Eq. (|A"5|)1. The lattice central potential can be fitted to the ansatz [including 
g of Eq. (fUD as a fit parameter] , 

VcoUR)-V c = KR-^ + ^ 2 -± (51) 

with self energy V c , string tension K and Coulomb coefficient e. The 1/R 2 and 1/R 3 func- 
tional form of the corrections corrections that account for the running of the coupling is not 
meant to be physical but has just been introduced to effectively parametrize the data within 
the given range of .R-values. For the j3 = 2.74 data we have set d — 0, such that we have 
5 fit parameters in this case while we allow for 6 parameters at (3 = 2.96. The resulting 
parameter values are displayed in Table |J. For technical reasons (link integration), only 
potential values for R > \/2 have been obtained, such that the fits do not include R = 1. 

In Fig. |, the potentials V COTT — V c from both /3-values are displayed in units of the 
string tension, K, as extracted from the fits, together with a fit curve that corresponds to 
the (averaged) values of fit parameters e = 0.257, / = af — 0.0054/a/k and d = a 2 d = 
0.00011/k. As can be seen, the two data sets scale nicely onto each other which means that 
the result appliesg to the continuum. Violations of rotational invariance are removed by our 
correction method, even at very small values of R, and the data is well described by the 
parametrization over the whole r range. 



B. Spin dependent Potentials 

Our results on the first spin-orbit potential V[ are displayed in Fig. [14] (in units of the 
string tension). The two data sets show approximate scaling behavior. In addition to a 
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constant long range contribution —K as expected from Eq. fll6|), we find an attractive short 
range contribution that can be fitted to a Coulomb-like ansatz, 

Vl(r) = -^ 2 -K. (52) 

For these one parameter fits we have constrained the constant long range part to the value of 
the string tension, as obtained from the central potential. We find the values h = 0.0474(58) 
and h = 0.0439(23) at (3 = 2.74 and f3 = 2.96 for this dimensionless parameter, respectively. 
Averaging these two numbers yields h = 0.0458(25). This additional term comes out unex- 
pectedQ and amounts to about one fifth of the e/r contribution to the central potential. 

Taking into account the running coupling improved effective parametrization of the cen- 
tral potential [Eq. (0)] and the Gromes relation, we expect 

V?M~-| + ^ <»> 

as opposed to the tree level expectation Eq. flT2]) . The Coulomb coupling has to be reduced 
by the amount that goes into V{. 

If we assume V{ to be generated exclusively from scalar-like exchanges and neglect the 
possibility of pseudoscalar contributions to V3, Eq. (§) leads us to modify Eq. ( |T3"D and 

expect, 

W = 5<i^-5f + ™. ,54) 



In Figs. [T| and 0, the spin-orbit potential V 2 and the spin-spin potential V3 are dis- 



played, together with the expectations Eqs. (B^) and (|54]), respectively. In both cases, we 



observe reasonable agreement between data and expectations. For small values of R, scaling 
violations between the two data sets from the different /3-values are evident as well as (in 
case of V3) some deviations from the expectation. As we will see in Section [IV C| , these dif- 
ferences can be explained as finite a effects and understood in terms of lattice perturbation 
theory. 

In Fig. the spin-spin potential V4 is displayed in lattice units for the two f3- values. An 
oscillatory behavior is observed which is similar to that of the lattice 5-function, expected 
on the tree level [Eq. ( |A20 ) of the Appendix]. Moreover, the two data sets nearly coincide 



with each other, in distinct violation of scaling. Higher order corrections to the 5-function 
which might scale with an appropriate dimension should account for the differences between 
the two data sets at small R. 



7 However, from the results of Ref. [20] and the Gromes relation, we can derive the following 
one loop connection between results on V{(r;fii), obtained at lattice spacings a\ = ir/fii and 
a 2 = 7r//x 2 @: V{{r;yi 2 ) = V[(r;ni) - (l - (a(/x 2 )/a(/ii)) 9 / 25 ) V 2 \r; W ), which means (i) that 
such a contribution must exist and (ii) that its relative weight will increase with decreasing lattice 
spacing. However, the actual magnitude of this admixture is still surprising as the prefactor of the 
admixture from V 2 is as small as about 0.03 under a scale change by a factor two. 
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C. Finite a aspects 



In Figs. 0-|2D|, we focus on the small R behavior of the sd potentials V^', V3 and V4. We 
restrict ourselves to display the (3 = 2.74 results only which are in qualitative agreement with 
those obtained at (3 = 2.96. In addition to the data (with error bars), the tree level pertur- 
bative expressions from the Appendix are displayed [squares, Eq. ( |A31| ) for K,', Eqs. ( |A23| ) 



and ( A24|) for V3 and Eq. ( A20|) for V4). The normalization constants c = Cpa s have been 



obtained from fits to the first seven data points. All three sd potentials are qualitatively 
described by these one parameter fits and deviations of the data from a continuous curve 
can be understood in terms of this lattice expectation. 

The fit parameters are displayed in Table |T|. From the analysis of the central potential, 
we would expect values c = e — h « 0.21 while the tree level lattice expectations [with lattice 
coupling a s — 1/ (7T/3)] are c = 0.087 and c = 0.081 for (3 = 2.74 and (3 = 2.96, respectively. 
In agreement with the perturbative expectation, all Cj come out to decrease with increasing 
(3. We find c = e — h as determined from Vq and V[ to be about 2.5 times larger than 
the naive tree level value while this factor reduces to 1.8 in case of and 1.3 for V3 and 
V4 as these potentials are dominated by higher momentum gluon exchanges and thus more 
perturbative. In order to investigate if the remaining differences between data points and 
tree level expectation (with renormalized coupling q as fit parameters) can be explained 
in terms of higher order perturbative corrections, we attempt to model running coupling 
effects. 

The only additional diagrams that contribute to Vq at 0{g i ) on the lattice as well as in 
the continuum are one loop corrections to the gluon self energy. The renormalization of the 
coupling, emanating from these diagrams, has been computed on the lattice for on-axis sep- 



arations of the sources ||33| , |34| . One can account for this correction by building in a running 
coupling constant a(q) into the gluon propagator in momentum space [Eq. ([A~5|)1. Instead 
of attempting to compute the correct lattice sum, we model this effect by the corresponding 
continuum expression, 

1/6 b 2 V 1 
a(t) = — l + -lnt + -lnt (55) 



Anb t \ t t 2 



with 



f 1 (?\ h 11N h 34iV h 

t = ln UV ' ° = ^ ' 6l = 306W ' W (56) 

where we replace q 2 by its lattice counterpart q 2 = 4J2i sin 2 (gj/2) to account for the periodic 
boundary conditions. The case 6 = corresponds to the one loop result. The difference to 
the correct lattice expression of Ref. |M] is small. A is a QCD scale parameter that can be 
related to the usual schemes via perturbation theory. 

To 0(g 6 ), apart from a renormalization of the gluon propagator, additional exchange 
contributions appear that can be resummed into a single running coupling by renormalization 
group arguments. These arguments do not apply to the lattice where rotational invariance 
is broken, such that, in addition to its absolute value, the direction of q has to be taken into 
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account. Bearing this in mind, we will nonetheless attempt to model two loop effects by the 
continuum two loop running coupling of Eq. (|55D . 

In case of the sd potentials V^-V^ not only the gluon self energy contributes to 0(g 4 ) 
but also exchange diagrams between the ears, incorporating a three gluon vertex. In the 
continuum these can be resummed into an effective running coupling. Due to this resum- 
mation, the scale parameters A, (for V^) can differ from each other. Again, we attempt 
to model this effect by plugging the continuum running coupling at scales q into the lattice 
tree level expressions. 

In case of the sd potentials, the infrared region is suppressed by powers of q [Eqs. ( |A38p - 
( |A40| )1, such that the form of the propagator in the non-perturbative domain has little 
effect. To remove the unphysical pole at q = A, an infrared protection can be build in into 
the propagator by substituting t by t c = ln(g 2 /A 2 + c 2 ) with a constant c. The smallest 
momentum on the finite lattice is q = it/ (aL a ). We choose c 2 = max(0, e — tt 2 / (aL a A) 2 ) (e 
is the Euler constant) to guarantee t > 1 with c 2 being negligible at large momenta q « 1/a. 

We do not attempt to fit the central potential to its perturbative expectation because 
this would require us to put a l/g 4 -like piece into the propagator by hand to generate the 
string tension. Also, as we found out, a large part of the short range structure is of scalar 
nature and would have to be modeled separately. A running coupling analysis of the central 
force in position space can be found in Refs. [ 15 ,[16 ] . 



Fits of the one loop as well as the two loop running coupling improved expectations (with 
one free parameter A) to the first 4-8 data points of each potential have been performed. 
The results of the two loop fits to 7 data points are included into Figs. [3~8| p0| (circles, dotted 
lines). The running coupling significantly improves agreement with the data. Differences 
between one and two loop results are small. The A parameters remain stable against the 
variation of the fit range within 10 %. 

In Table |TTT| , results on the A parameters from our one and two loop fits are collected. 
We do not include any errors since the fits are only thought to qualitatively describe the 
data with reduced x 2 - va lues of typically 10-100. We observe approximate scaling between 
the two sets of A parameters at (3 = 2.74 and j3 = 2.96 within 15 %. However, the two loop 
values are larger by about a factor two than the corresponding one loop values. 



V. CONCLUSIONS AND OUTLOOK 

We have devised methods to determine spin dependent interquark forces to high pre- 
cision. From our high statistics lattice simulation in SU(2) gauge theory, we find reliable 
renormalized potentials with good scaling behavior. There is clear evidence for a short range 
scalar exchange contribution in the long range spin-orbit potential at the level of 20 % of 
the Coulomb part of the central potential. The other sd potentials are found to be short 
ranged and are well understood by means of perturbation theory. 

An extension of the present investigations to the case of interest, SU(3) gauge theory, 
and inclusion of velocity dependent corrections will be presented in II ||18j| . As a further 
step, predictions from various models of QCD interactions, such as dual QCD |55j can be 



tested against lattice results on the potentials. Also spectra, wave functions and decay 
constants for arbitrary (sufficiently large) values of the quark masses can be computed just 



18 



by solving a simple differential equation f36| . Subsequently, these results can be confronted 
with experiment or compared to results from direct lattice NRQCD predictions || as a first 
principles check of the viability of the instantaneous approximation of the potential picture. 
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APPENDIX A: WEAK COUPLING EXPANSION OF LATTICE POTENTIALS 

In this appendix, we will elaborate the tree level expectations for the sd potentials from 
a weak coupling expansion of SU(N) gauge theory on the lattice. Some of these expressions 
have already been derived in Ref. ]13[ for on-axis source separations. To illustrate the 
method, we start with the central potential before considering sd corrections. We include 
the corresponding continuum expressions for completeness. 



1. Central potential 

A weak coupling expansion of the Wilson loop yields to lowest order in g 2 , 

(W(R,T)) = ex V (-Cpg 2 ~ ^WG(m-n)y,(n) (Al) 



T-l 



exp C F g 2 T £ [G(R, t) - G(0, t)] (T » R) (A2) 



t=o 



where J M (m) = ±1 if (m, fi) G dW and J^(m) = elsewhere. Only terms, extensive in T, 



have been kept in Eq. (|A"2|). For SU(2), the color factor Cp = 2N becomes Cp = 3/4. 
G{n) denotes the lattice gluon propagator in position space, 

G H = T^-E^ , ^ = 2sin(|) (A3) 



with 

Qi = j^rrii , mi = — - + 1, . . . , — , 

27r L T L T i . 

q 4 = — m 4 , m 4 = — - + 1, . . . , — . (A4) 

Note, that we have neglected the zero momentum contribution to W which is suppressed by 
a factor RT / (L^L T ). 
With 
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G*(R) = £G(B,f) = J_E V 

(J q^O ^ 



D iqR 



(A5) 



and Vq(R) = — lim^oo In (VK(R, T)) /T, one obtains 



1 



K)(R) = -C F g 2 (G L (R) - G L (0)) — >V C - C F a s - (R - oo) 



(A6) 



for the central potential with a s = g 2 / (47r). V c = C F g 2 GL(0) denotes a self energy constant. 

2. Spin-spin potentials 

To compute V{(H) - V^(R) from Eqs. ( p~7| ) - (p0|), only gluon exchanges between the two 
color field insertions (ears) have to be considered to order g 2 . Starting from 



(P^) = exp(-C F g 2 X_ 



12 



with 



n 



(A7) 



(A8) 



where = ±1 (the sign depends on the orientation of the link) if (n, /i) G dPi and 

J^(n) = elsewhere, we obtain to order g 2 [Eq. (|22|)1, 



<*> 


w (p 1 -p!)(p 2 


- ph_ 


)(W) 




w (Pi + Pi) 


){v 


W {P 2 + P ] 2 ) 


) 



\({Pi-pI){P2-pI: 



(A9) 



sinh (C F (? 2 X 12 ) = C F 2 X 12 , 



(A10) 



with Fj 



2i 1 J r j 



Note, that the self energy contributions from interactions within 
one ear are not small, though of order g , but a 60 % contribution at our /3-values. However, 
such contributions are cancelled by the denominator of Eq. (|A9|). Within the potentials V^' 
- V4, all order g 6 self interactions are cancelled as well. 

Let us consider the correlation function between two magnetic ears at positions (0, 0) 
and (R, t) where we choose the first ear to be within the i — j plane and the latter within 
the i — k plane with i 7^ j, i 7^ k. In this case, we obtain 

V 12 = 2j o °°dt ((F&Vw = -C F g 2 (A^A^ + <^AHa| +) ) G l (R) (All) 



where 
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A l (+) /(n) = /(n + i)-/(n) , ^/(n) = /(n) - /(n - i) (A12) 
are forward/backward differences. Note, that 

("OO 

G L (R) =2 / dtG(R,t). (A13) 



o 



Interactions between averages of four adjacent plaquettes centered around lattice points 
and R, which we use for the magnetic ears [Eq. fl24|)], can be derived from Eq. ( All ) by 
averaging over the 16 possible combinations of single plaquette ears, 

V£ = -C F g> (A<">A« + S jk ^M +) ) S^sSr^CR) (A14) 

= -C F g 2 (A^-AfcSi + <y ifc A|Sj) Gl(R) (A15) 

with 

2 l (+ V(n) = ^(/(n) + /(n + i)), (A16) 

S l H /(n) = i(/(n) + /(n-S)), (A17) 

Si/(n) = i (2/(n) + /(n + i) + /(n - i)) , (A18) 

Ai/(n) = i(/(n + S)-/(n-S)). (A19) 

Note, that Hj = l(S< + > + 2«) = sf^H, E,Ej = 1(2, + -,), A, = i(A, (+) + A«) = 

A^ + 'El| "* = A- Since all H and A are linear combinations of translations, they 

commute with each other. 

In case of V4 the two magnetic ears are parallel, i.e. j = k, such that a sum over the 
three possible i,j combinations [Eq. fl20|)1 yields, 

t/ 4 (R) = -2C F g 2 ]T A^sf'G^R) = -ICrfl&E G L (R) (A20) 



i=l 



with 



, 2 = 11:2*, (A21) 



Af)=A«A(-) , A( 2 ) = £AS 2) . (A22) 

Note, that A? = A^Hj. 

Correlators between Bi and _B m are required for computation of V 3 . Since the direction 
of B is orthogonal to the plane of the corresponding plaquettes, we find the relations, k = I, 
j = m and i 7^ j, i 7^ k for I 7^ m. In case are cyclic, (m,i,k) are automatically 

anticyclic, such that we obtain an overall minus sign. We find for j ^ k from Eq. (|l9l)Fl, 



i Obviously, this expression is only useful for off- axis separations where R~ 7^ and 7^ 0. 
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R 2 

For I = m = i we obtain: 



V 3 (R) = Cpg'AjA^GLiR). (A23) 



3R U V 3 = C F g 2 (2A 2 E ( i ±) - A 2 (2~ k - - A^(2S i - G L (R) 

= C F , 2 i (5Af E - A^H - Af% - AfE,) G L (R) (A24) 



with j ^ i, k ^ i, j ^ k, R^ = RiRj/R 2 — (5^/3. 

It is easy to see that the above expressions amount to 



and 



or 



VAr) = -2C F a s V 2 - = 8nC F a s 5 3 (r) (A25) 
r 

V 3 (r) = Cpas—d^- = 3C F a s \ (A26) 
r^j r r 6 

V 3 (r) = C F (R u y l1 - (2d 2 -d 2 - d 2 ) l - = C F a s (R u y l ^t^ = 3C F a s ^ (A27) 
in the continuum limit. 



3. Spin-orbit potentials 

For computation of V[ and V% one has to take into account correlators between plaquettes 
in the i—j plane and the i — 4 plane. To lowest order in g only exchanges between the links 
oriented in % direction have to be taken into account. With 

roo 

2 dtt [GCR, t) — GCR, t + 1)1 = G L (R), (A28) 
Jo 

we obtain 

roo „ _ 1 , . 

W 12 = / dtt((F 1 F 2 )) w = -C F g 2 -A^G L (R) (A29) 
Jo 2 

for the integrated correlation function. Averaging over the relevant plaquette combinations 
finally yields, 

W& = -CV^A^G L (R). (A30) 

Thus, V[ vanishes to order g 2 while the leading order expression for the second spin-orbit 
potential [Eq. flT8|)] is, 

^(R) = - jr C^AjE^ G L (R) . (A31) 

JXj 

In the continuum this amounts to, 

V&r) = -Cpas-dj- = C F a s \. (A32) 
Tj r r 
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4. Tree level relations 



Tree level relations between the lattice potentials that are analogous to Eqs. (0) and (Bp 
can be derived: 

T/ 2 '(R) = I-A.-Sf >Vfc(R) + i?(R), (A33) 
y 4 (R) = 2Af } H t> (R) = 2 £ A, (f^'(R)) • (A34) 

From Eq. ( |A33| ) one might attempt to generalize the Gromes relation. Let us assume for 
the moment that a linear difference operator exists, such that 

V£(R) - V((R) = £ c(n)K,(R + n) (A35) 

n 

with constants c(n). Both sides of the above equation can be expanded in orders of g 2 . At 
order g 2 we find, 

5>(n)G L (R + n) = ^A,sf } G L (R). (A36) 

The factor R/Rj illustrates that c(n) has to depend on R, in contradiction to the ansatz, 
i.e. nonlinear corrections have to be included. Also, in our numerical studies we find the tree 
level relation Eq. ( |A33| ) to be substantially violated at small R. Of course, the continuum 
Gromes relation as well as the above lattice version are retrieved at large R. 



5. Continuum results 

For continuum potentials one obtains the following tree level expressions, 



Vn(r) 



f dq 3 e iqr 2a s r°° sin qr a s 

-C F a s / = ~C F / dq = -C F —, (A37) 

J 2ir z q z it Jo qr r 

Vi{r) = -iC F a s I ^^> r = -C F ^ (°° dqq 2 r j^qr) = C F ^, (A38) 
J 27T Z q z r 7C Jo r z 

V 3 (r) = -C F a s I ^-^e^ = -C F ^ rdqq 2 j 2 (qr) = 3C F ^, (A39) 
J 2tt z q l r A it Jo r 6 

VAr) = C F a s [ ^e iqr = C F ^ f°d gg 2 ^^ = 87iC F a s 5 3 (r). (A40) 
J 2n 2 7r jo qr 

A linear confining contribution can be introduced by adding a 1/g 4 term to Vq in momentum 
space. The integrals for the sd potentials are suppressed at low q like q 2 or q 3 , such that 
we naively expect perturbation theory to be more reliable in this case than for the ground 
state potential. Also, finite size effects are expected to be smaller. V[ vanishes at tree level. 
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TABLES 



TABLE I. Simulation parameters. The physical scale has been obtained from 


/ii = 440 MeV. 


P 


V = L\L T 


a/fm 




L a /fm 




^■conf 


2.74 


16 4 


0.041 




0.65 




1290 


2.74 


32 4 


0.041 




1.30 




200 


2.96 


32 4 


0.022 




0.70 




204 


2.96 


48 4 


0.022 




1.05 




67 


TABLE II. Fit parameters 


to the central potential [Eqs. (|49|)- 


mi 


Parameter 


P = 2.74 




P = 2.96 






Average value 


K 


0.00785(13) 




0.00246 (8) 






— 


e 


0.2611(35) 




0.2507(41) 






0.2567(74) 


/ 


0.0598(33) 




0.1114(76) 






n nn^zii f^s^ u/~w 

u.uuo^i yoo j / v 


d 






0.0437(50) 






0.000108(12)/^ 


Vc 


0.4898(12) 




0.4334 (9) 








9 


0.1707(12) 




0.1409 (9) 








TABLE III. 


The tree level constants c 


: = Ci?a s and A parameters from the 


running coupling 


analysis of V 2 ', V3 and V4. 












tree level 












c 4 


P = 2.74 


0.17 




0.12 






0.13 


P = 2.96 


0.14 




0.10 






0.11 


running a 






K 




A 4 / v / k 


1-loop P = 2.74 


0.29 




0.12 






0.08 


1-loop P = 2.96 


0.25 




0.10 






0.07 


2-loop P = 2.74 


0.67 




0.32 






0.20 


2-loop P = 2.96 


0.61 




0.29 






0.18 
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FIGURES 
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FIG. 1. Nonperturbative renormalization prescription for eared Wilson loops. 
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FIG. 2. Corrected central potential V COTT — V c in units of the string tension at (5 = 2.74 and 
(3 = 2.96. The fit curve corresponds to the parametrization Eq. ( |5T|) with parameter values as in 
the last column of Table |D|. 
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FIG. 3. Ground state overlaps as a function of the quark separation R (in lattice units) at 
P = 2.74. 
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FIG. 5. The correlation function (^(i) for R = 2. The solid curve denotes our interpolation 
while the dashed curve corresponds to a three parameter fit to the data. 
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FIG. 6. The correlation function C2(t) for R = 4. The solid curve denotes our interpolation 
while the dashed-dotted curve corresponds to a three parameter fit to the data. 
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FIG. 7. The correlation function C^{t) for R = 2. The solid curve denotes our interpolation 
while the dashed curve corresponds to a four parameter fit to the data. 
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FIG. 8. Test of the Gromes relation Eq. (||). The combination V{ — V[ is compared to the 
central force as obtained from the parametrization Eq. (|5l|). 
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FIG. 9. Comparison of the central potential Vq(R) at (3 = 2.96 between a 32 4 and a 48 4 lattice 
(in lattice units). 
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FIG. 10. Comparison between —V{(R) at (3 = 2.74 obtained on a 16 4 lattice and a 32 4 lattice 
(in lattice units). 
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FIG. 11. Comparison between V^R) at j3 = 2.74 obtained on a 16 4 lattice and a 32 4 lattice (in 
lattice units). 
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FIG. 12. Comparison between V^{R) at (3 = 2.74 obtained on a 16 4 lattice and a 32 4 lattice (in 
lattice units). 
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FIG. 13. Comparison between V^{R) at (3 = 2.74 obtained on a 16 4 lattice and a 32 4 lattice (in 
lattice units). 




FIG. 14. The spin-orbit potential V{, together with a fit curve of the form — V/(r) 
(with h = 0.046) in units of the string tension. 
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FIG. 15. The spin-orbit potential V£ in comparison to the continuum expectation from Eq. (|53|). 
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FIG. 16. The spin-spin potential V3 in comparison to the continuum expectation from Eq. (|54[). 
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FIG. 17. The spin-spin potential V4 for the two (5- values in lattice units. 
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FIG. 18. Comparison of the lattice potential V^' (points with error bars) to tree level (open 
squares) and two loop running coupling improved (open circles) lattice perturbation theory. 
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FIG. 19. Same as Fig. Ill for % 




FIG. 20. Same as Fig. 18 for V4. 
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